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We study the rapidity spectra in ultra-relativistic heavy ion collisions in the framework of the 
Landau hydrodynamical model. We find that thermal smearing effects improve the agreement with 
experimental results on pion rapidity spectra. We describe a simple model of the hadronization and 
discuss its consequences regarding the pion multiplicity and the increasing entropy condition. 
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PACS numbers: 25.75.Nq, 24.10.Nz 

I. INTRODUCTION 

The Landau hydrodynamical model [H, is a simple 
approximate solution of the hydrodynamical equations 
describing the expansion of a thin disk of static gas. The 
original idea was to use hydrodynamics to describe colli- 
sions of very high energy hadrons, when a large number 
of new particles are created. Today the model is also 
frequently applied to the expansion of the Quark Gluon 
Plasma (QGP) formed in ultra-relativistic heavy ion col- 
lisions. It provides a good description of the shape of the 
rapidity distribution of pions [3j . 



II. REVIEW OF THE LANDAU MODEL 

The initial state of the Landau hydrodynamical model 
is a thin static disk of thickness A and diameter a. It is 
assumed to be the approximation of the overlap of two 
highly Lorentz contracted nuclei, and the two geometrical 
sizes are assumed to be related by 



A = a/ 7 , 



(1) 



with 7 the Lorentz factor corresponding to the velocity 
of the colliding nuclei in the center-of-mass frame, 
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where Snn is the invariant collision energy per colliding 
nucleons of mass m n . 

The hydrodynamical evolution of the system is divided 
in two parts: the initial longitudinal expansion (dur- 
ing which transverse velocities and displacements are ne- 
glected), and the subsequent "conic flight", where trans- 
verse velocities appear. 

Assuming a simple equation of state (EOS) of the form 



P = e/3, 



(3) 



lj. This solution is summarized (using the notation of 
3]) as follows. The energy-density field is 



where y± are logarithmic light-cone coordinates 

t±z 



y± = in ■ 



A 



(4) 



(5) 



and eo is the initial energy-density of the disk. The ra- 
pidity field is 

y(t,z) = = \{v+-y-)- ( 6 ) 

In the longitudinal expansion components of the flow 
four- velocity are expressed in terms of the flow rapidity 

as 

= (coshy, 0, 0, sinhy) = coshy(l, 0, 0, v z ) (7) 
which can be inverted to give 
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(8) 



Looking at ^ one observes that the flow rapidity and the 
space-time rapidity coordinate coincide (the first expres- 
sion in (|6]) is exactly the definition of space-time rapid- 
ity), similarly to the case of the Bjorken model. Because 
of the coincidence of flow rapidity and space-time rapid- 
ity the space-time coordinates of a fluid element and the 
components of its four- velocity are related by 



(9) 



where we have introduced the proper time coordinate 

T = Vt 2 " Z 2 . 

In Refs. P, IH an estimate of the transverse displace- 
ment of the system is given as 



an approximate solution of the 1+1 dimensional problem 
of the longitudinal expansion phase is given by Landau 
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The transition to the second phase of "conic flight" is 
assumed to happen when this transverse displacement 
reaches the transverse size of the system, a. Using the 
relations Eq. ([5]) the condition of transition to conic flight 
can be written as 



ttr 



2a. 



(11) 



This defines a hypersurface in space-time, which we call 
the transition hypersurface. 

Although 0,1] give a sketch of the hydrodynamical 
evolution of the system in the second phase of conic 
flight, in the calculation of observables they assume that 
freeze-out happens immediately, i.e. not only the flight 
angle 6 of fluid elements is assumed to be frozen on the 
transition hypersurface defined by Eq. ([TT]) . but also the 
fluid is replaced by an ensemble of non-interacting par- 
ticles traveling with the same constant velocity. Thus, 
the transition hypersurface is also the freeze-out (FO) 
hypersurface (and the subscripts TR and FO can be in- 
terchanged). The velocity of particles is assumed to co- 
incide with the velocity of the fluid element at the FO 
hypersurface, which means that no smearing effect com- 
ing from the thermal distribution of particles in the fluid 
is taken into account. 

Equation ([TT]) gives a relation between the time of 
freeze-out and the position of the fluid element at freeze- 
out of the form 



t FO (z) = t TR (z) = \/4a 2 + z 2 . 



(12) 



This means that at the freeze-out the system can be de- 
scribed in terms of functions of one variable, say zpo- 
(Because of azimuthal symmetry everything is indepen- 
dent of the azimuth angle (f>, furthermore, everything is 
assumed to be independent of the radial coordinate r for 
r < a, and to vanish for r > a). Furthermore, there is a 
one-to-one correspondence between z and the fluid rapid- 
ity y on the FO hypersurface of the form (see Eqs. ([7]), 
© and ([IT])) 



zpo = 2asinhy, 



(13) 



so instead of zfo we can use y as the independent vari- 
able. Using Eq. (|T2|) the FO time can be expressed in 
terms of y as 



tpo = 2a coshy. 



(14) 



In order to get the rapidity distribution of particles we 
start from the expression of the total particle number 



N = / N»d<Ta, 



(15) 



with — nu^ the particle current (n is the invariant 
scalar particle density), and da^ the hypersurface ele- 
ment four-vector. Using cylindrical coordinates (r, <f>, z) 
the FO hypersurface is given by 



x^{r,(f),z) — (t-po{z),r coscj), rsincf), z) 



(16) 



and so 

da^ = e^pg-drx" d t px p d z x cr drd(f)dz = (1, 0, 0, — d z t-po(z)) rdrdtfidz. 

(17) 

From Eq. (TTT 



d z t FO (z) = tanhy, (18) 
and changing the integration variable z to y gives 

dz — 2a coshydy. (19) 
The <p and r integrations give a factor 



2tt pa/2 
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rdr 



(20) 



so finally 



do> ->• -y- (cosh y, 0, 0, - sinhy) dy. (21) 

We now assume that the transverse component of the 
velocity is negligible, so u° and u 3 are given by ([7]) at FO 
too. Therefore 



u^da^ -> — dy. 



(22) 



We assume that QGP is an ensemble of massless quarks 
and gluons, and obeys the Stcfan-Boltzmann EOS, e oc 
T 4 (which is consistent with ([3])). In that case 

n oc T 3 oc e 3/4 = e^ 4 exp {- (y + +y_ - y/y+yJ) } , 

(23) 

where we have used Eq. (UJ) in the last step. As described 
m @, from Eq. © it follows that on the FO surface 
U±\fo =yb±V, where y b = ln(2a/A) = ln(y/s NN /m N ) 
is the rapidity of the colliding nuclei in the CM frame. 
Therefore the exponent in (|2"3"|) can be written on the 
FO surface as (y+ + y_ - y/yyy^) |fo = 2y b - \Jvl~ V 2 - 
Using this in (|23[) we get the invariant particle number 
density on the FO hypersurface expressed as a function 
of y: 



n(y) ocexptJy 2 - y 2 



(24) 



The rapidity distribution of particles, dN /dy is obtained 
by differentiating (TT51) . Thus, making use of (|2"2"j) and 
([541 we get 



dN/dy ex exp j \Jy 2 b -y 2 ^ , (25) 
in agreement with 

III. THERMAL EFFECTS 



The starting point of the calculation of observables is 
the momentum distribution of emitted particles, which 
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develops at freeze-out. In the case of sudden FO hap- 
pening on a hypersurface the invariant momentum dis- 
tribution is given by the Cooper- Frye formula [IJ 



d 3 N 

= I f(x,p)p^da fl , 



(26) 



where f{x,p) is the phase-space distribution of particles 
after freeze-out. 

In hydrodynamics local thermal equilibrium is as- 
sumed, therefore the x dependence of / is encoded in the 
space-time dependence of the parameters of the thermal 
distributions, T(x) and u^(x). (The Stefan-Boltzmann 
EOS e oc T 4 restricts the considerations to baryonfree 
media, therefore the baryo-chemical potential is zero.) 
The thermal phase-space distribution has the general 
form 



f(x,p) = f 



c ( P^U^X) 



T(x) 



1 



(2tt) 3 



exp 



V T(x) 



(27) 

with £=1 for fermions, -1 for bosons and for classical 
particles obeying Boltzmann statistics, g is a degeneracy 
factor that is different for each particle species. 

In the case of the Landau model parameters of the 
fluid taken on the FO hypersurface depend on only one 
variable, which can be chosen the fluid rapidity y. The 
expression for T(y) can be obtained from the expres- 
sion of the energy density, Eq. assuming the Stefan- 
Boltzmann EOS e = <jT a . Following the steps leading to 
the expression ([25]) we get the result 



!/4 e- 2 W3 exp jl 



y'b 



V 



(28) 



The transverse velocity of fluid elements is assumed to 
be negligible at FO. Components of the four- velocity of 
the fluid are given (in accordance with ©) by 



u^(y) — (coshy, 0, 0, sinh y). 



(29) 



A comment concerning the validity of the model is ap- 
propriate here. Because of the relation u^u^ — 1 the in- 
clusion of a nonzero transverse velocity of the fluid would 
require the modification of Eq. ||7J) . But the relations J7]) 
are a key ingredient of the Landau model related to the 
separation of longitudinal and transverse expansion of 
the system. Their modification would destroy the sim- 
plicity of the model. 

An estimate of the transverse velocity of fluid elements 
can be obtained by differentiating (fTU|) . We get 



dr(t,y) 
Using Eq. (fl4|) this gives 

V r \ F Q = 



t 



2a cosh y 



coshy 



(30) 



(31) 



on the FO hypersurface. On the other hand Eq. (J7J gives 
v z = u 3 /u° = tanhy for the longitudinal velocity com- 
ponent, which implies 



vl = 1, 



(32) 



meaning that fluid elements would reach the speed if 
light at FO if one included a transverse velocity with 
the above approximations. This also indicates that the 
Landau model — although it gives a good description of 
rapidity spectra — is not a suitable model for describing 
e.g. the transverse momentum distribution of particles, 
at least not in its simplest form. 

The four-momentum of particles can be specified as 

p M = (m,T cosh y p ,pr cos 4> p ,pt sin <p p , vtlt sinh y p ) (33) 

in terms of their transverse momentum px, rapidity y p 
and azimuth angle <fi p , where tot = \J m 2 + p% is the 
transverse mass, with the mass of the particle, to. From 
Eqs. (J2SJ and IPH 



p^u^x) = TO T cosh(y - y p ). 



(34) 



In order to calculate the rapidity distribution of emit- 
ted particles one first has to express the invariant mo- 
mentum distribution, Eq. (|26|) in variables px, y P and 



E- 



d 3 N 



d 3 N 



dp 3 prdpxdypC 



(35) 



Integrating this, and substituting the right hand side of 
(|26p we get for the rapidity spectrum 

dN f d 3 N f f 

— = J dpTPTd4> p E-^ = J dp T pTd(j>p J f(x,p)p' 1 da lJ 

P .... ( 36 ) 
Here da^ is given by (|17p . The integrand is independent 

of the azimuth angles <p an d 4> P corresponding to the po- 
sition of the fluid element and the particle momentum, 
and the radial coordinate r (see Eqs. (gTJ), (® and 
This means that the effect of the da^ integration is the 
same as given by (|2"TT) and the integral over <fi p yields 
a factor of 2ir. The rapidity distribution is then (also 
making the change of integration variable px — > tot) 



dN 
dy P 



drriTm3p J dy cosh(y 



V P )f 



7 (P^u p {y) 



T(y) 



(37) 



IV. 



THE HADRONIZATION PHASE 
TRANSITION 



In accordance with the Stefan-Boltzmann EOS as- 
sumed in the Landau model the expanding quark-gluon 
plasma (QGP) is assumed to be an ensemble of massless 
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quarks, antiquaries and gluons. The EOS also requires 
that the QGP is baryonfree, which is a reasonable ap- 
proximation at RHIC or LHC energies, where the particle 
spectra are dominated by secondaries. 

On the other hand, the post freeze-out (PFO) hadronic 
medium is an ensemble of massive particles of different 
species. We assume that the phase transition from QGP 
to hadronic matter happens at the FO hypersurface (si- 
multaneous FO and hadronization). Physical parame- 
ters of the two media are related via the prescription of 
boundary conditions on the FO hypersurface describing 
energy and momentum conservation, 

[T^da v ] = 0, (38) 

where the meaning of the square bracket is [a] = o|pfo — 
a | qgp • The energy- momentum tensor has the form 

= (e + P)u^u u - Pg^ (39) 

in ideal hydrodynamics, da^ is unit four-vector normal 
to the FO hypersurface. In the Landau model this is 
given by 

da^ = da^j \J daPda p = (coshy, 0, 0, — sinhy) = 

(40) 

(see (f2"Tj) and ([29"]) ). The last equation reflects a gen- 
eral property of the Landau model, namely, that the 
four- velocity of the fluid is everywhere normal to hyper- 
surfaces characterized by a constant space-time rapid- 
ity r, in particular to our FO hypersurface. This fea- 
ture of the Landau model is the same as in the Bjorkcn 
model. Making use of relation (|40p the condition of 
energy-momentum conservation Eq. (|38[) becomes 

feu"] = 0. (41) 

The component of (|4lT) parallel to da^ tells us that the 
invariant energy density is constant across the FO hyper- 
surface, [e] = 0, while the disappearance of the compo- 
nent normal to da^ implies that the four-velocity of the 
fluid is unchanged during the phase transition. 

Although the invariant energy density remains the 
same at FO, the EOS changes from QGP of massless 
quarks and gluons to hadronic matter of massive con- 
stituents. This causes a change of temperature even if 
the FO normal, <i<7 M , and the local flow velocity, u^, co- 
incide, so that the flow does not change. 

In addition to energy-momentum conservation one also 
has to check the validity of the increasing entropy condi- 
tion, which yields a boundary condition of the form 

[s«1 > 0, (42) 

where s is the invariant entropy density. 

In order to reduce the numerical complexity of the sys- 
tem we assume that the massive particles in the hadronic 
phase obey the Juttner (relativistic Boltzmann) distri- 
bution. This approximation is justified at the typical 



phase-transition temperatures of T sa 170 MeV. Then 
the pressure of the medium is given by 

* — ' 7T 

a 

where m a and g a are the mass and degeneracy factor of 
particle type a, and Ki is a modified Bessel function of 
the second kind. The energy-density is 

where we used the notation 

K la =K l (m a /T). (45) 
Finally, the entropy density of the system is 

s ( r ) = E^=E^ m «K^ + 4T ^)- (46) 

ol a 

The temperature of the medium at a given position on 
the FO hypersurface is obtained by solving Eq. (|4"4")) nu- 
merically. Then the phase-space distribution function, 
Eq. (|27|) with £ = 0, is known on the FO hypersurface, 
and observables can be calculated via integration - e.g. 
the rapidity distribution as given by Eq. (|3"6")l . 

Now we have to give the explicit form of the EOS in the 
hadronic phase by specifying the particle types present 
in the medium and fixing their properties (mass and de- 
generacy). At ultra- relativistic energies most of the pro- 
duced particles are pions. Based on this we can model 
the hadronic phase by a pion gas with m w = 0.139 GeV 
and g-x = 3. One expects that this model will overesti- 
mate the pion yield because in reality some of the energy 
is distributed among other hadron species. However, this 
model can be used as a reasonable first approximation. 

Below we briefly discuss two other possible models for 
the post FO medium. 

At RHIC the constituent quark number scaling of flow 
observables has been found Q- This observation is con- 
sistent with simple versions of the coalescence (or recom- 
bination) model of hadronization (see e.g. @). This in 
turn suggests that the collective motion of the medium 
is developed already at the constituent quark level, and 
is basically unchanged during the process of recombina- 
tion. In practice this means that the medium can be 
modeled by an ensemble of constituent quarks and an- 
tiquarks when one studies collective observables such as 
momentum spectra of abundant particles. 

In this approach pions need a special treatment be- 
cause - due to their small mass - they can not be de- 
scribed by the recombination of a constituent quark and 
antiquarks of mass ~0.3 GeV. Indeed, pions are special in 
the sense that they are the Goldstone bosons associated 
with the breaking of chiral symmetry, which is responsi- 
ble for the generation of hadron masses. 

Based on the above we can model the post FO hadronic 
medium by an ensemble containing pions in addition 
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to constituent u, d and s quarks and their antiquarks. 
We will assume the values m U) <j = 0.31 GeV and m s — 
0.5 GeV for the constituent quark (CQ) masses, and use 
the degeneracy factors g u = ga = g s = 12 where we 
have taken into account the number of colors N c = 3 the 
spin degeneracy factor of 2, and another factor of 2 for 
the inclusion of antiquarks. We assume that the explicit 
inclusion of pions does not affect the degeneracy of the 
quarks. This is justified e.g. if pions are formed earlier 
than other hadrons. 

As another possibility, the medium can be modeled 
by a hadron gas containing low mass hadron states, like 
pion, nucleon, kaon, r), p, and lo mesons, and A particles. 



V. PION RAPIDITY DISTRIBUTION 

The BRAHMS collaboration at RHIC has measured 
the pion rapidity distribution in central Au+Au collisions 
at ^/s^n = 200 GeV 0,9. In order to obtain the same 
spectrum from Landau hydrodynamics first one has to 
fix the parameters of the initial condition. The radius of 
the gold nucleus is taken to be r = 6.98 fm, the Lorentz 
contraction factor is 7 = \/sj\nv/2mjv = 106.6. These fix 
the values of a = 2r and A = a/7. The initial energy 
density of the system is estimated as e = E/V, where 
V is the volume of the Lorentz contracted gold nucleus, 



V = £r»/ 7 . 



(47) 



Note that the volume of the initial disk of the Landau 
model is 



V L =r 2 nA = 2r 3 /7 = 3V/2, 



(48) 



which means that the total energy content in the initial 
state of the Landau model is a factor of 3/2 larger than 
the energy of the reaction in reality. 

A comparison of the pion rapidity spectra obtained 
from the model and the experimental results can be seen 
on Fig. [TJ Three theoretical curves are shown. The first 
one (labeled "no T") corresponds to Eq. ([25]) . This is 
identical to the result of Q and was obtained by neglect- 
ing thermal effects. The second one ("QGP") shows the 
result of the calculation taking into account the thermal 
smearing in the QGP. This curve is the result of the in- 
tegral (f3"Tf where / is the Fermi distribution of massless 
quarks. Finally, the third one ("pion, PFO") is calcu- 
lated from the thermal pion distribution function valid 
after the phase transition, assuming a pion gas in the 
post FO side. In the case of the second and third curves 
the integral (|37p was calculated numerically. 

Looking at the shape of the obtained rapidity spec- 
tra we see that thermal effects improve the precision of 
coincidence of theoretical and experimental results. 

The author of @ discusses a few possible corrections 
to the Landau model. One of them is related to the 
fact that the initial compression of system depends on 
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FIG. 1: (Color online) Rapidity distribution of charged pions 
emitted in an Au+Au collision at ^Jsnn = 200 GeV. The 
experimental data points are from [3]. See the text for the 
explanation of the three theoretical curves. 



the EOS, therefore the thickness of the initial state can 
deviate from the value obtained by taking into account 
solely the Lorentz contraction of the nuclei. The other 
possible correction comes from the observation that an 
arbitrary multiplicative factor can be inserted in the con- 
dition of transition to conic flight, Eq. (fTTj) . since this 
equation expresses only a rough equality of the trans- 
verse displacement and the transverse system size. Both 
corrections lead to the modification of the rapidity dis- 
tribution, Eq. (|25|) . of the form 



dN/dy cx exp { VQ/b + C) 2 -y 2 } , 



(49) 



where the unknown quantity £ can be used to improve 
the agreement of the model and experiment by fitting to 
the experimental data. Without questioning the validity 
of these arguments we note here that some of the discrep- 
ancies of the Landau model and the experimental data 
can be explained by the inclusion of thermal effects. 

In the case of all three curves on Fig. HQ an overall nor- 
malization factor is used following [lj, |3j as a free pa- 
rameter and is fitted to the experimental results. This 
normalization factor, A/", is defined by the equation 



dN 
dy P 



1 dN 

normalized ™ ^Vp 



(50) 



primary 



where (dN '/ 'dy p ) primary is given by Eq. fl37j|. In the case 
of the first two curves ("no T" and "QGP") AT has no 
meaning, since the models do not differentiate between 
different particle species, and give no prediction for the 
total pion yield. In the third case ("pion, PFO") on 
the other hand the results are obtained from the pion 
phase-space distribution function, and the value of the 
normalization can be used to check the validity of the 
model concerning pion multiplicity. 
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Using the simple pion gas model for the post FO 
medium a normalization factor of JV = 9.29 was ob- 
tained, the primary model overestimates the data by this 
factor. [l(| The total four-momentum of the system can 
be calculated as 

PL = J T^da v , (51) 

where the integration is over an arbitrary spacelike hyper- 
surface, and the energy-momentum tensor T^" is given 
by Eq. ((35} with e, P and as in ©, fl3J, and ©, re- 
spectively. Because of the approximations used in the 
model the integral Eq. (|5ip is not independent of the 
choice of the hypersurface. Specifically, on a surface with 
constant proper time, r the Oth component of (|5ip . the 
total energy, has the form 

Ptot = E tot = —t / ecoshydy, (52) 

where in the integration we have used steps similar 
to those leading to Eq. (f2T|) . Calculating the integral 
Eq. (1521) numerically on the FO hypersurface character- 
ized by t — Tpo = 2a (see Eq. (fTTI) ) for an Au+Au 
collision at y / ijv r /v = 200 GeV, we get a factor of 5.16 
larger value for the total energy in the primary model 
than the energy of the reaction, A^/sjvjv. An overesti- 
mation of the pion yield with the same factor is a natural 
consequence. 

In the experimental analysis of @ the 5% most cen- 
tral events have been used, while the Landau model pic- 
tures an exactly central collision. Based on a simple geo- 
metrical model describing the colliding nuclei as spheres 
with uniformly distributed nucleons the average number 
of participants in this experimental sample is 88% of the 
number of participants in a central collision. Assuming 
that the pion multiplicity scales with the number of par- 
ticipants this means that a model valid for central colli- 
sions should give a factor of 1 /0.88 = 1.14 larger result for 
the pion multiplicity. Thus, centrality selection and the 
violation of energy conservation by the Landau hydro- 
dynamical model together explain a normalization factor 
of 5.16 • 1.14 = 5.88 in contrast to the value M = 9.29 
found in the fit to the experimental pion spectra. The 
remaining discrepancy can be attributed to the fact that 
the simple pion gas model of the past FO medium can 
not account for the energy carried by other hadrons. We 
can conclude that (5.88/9.29) • 100% = 63% of the energy 
is carried by pions in the final state of the collision. 

In Q ratios of the full phase-space extrapolated 
kaon/pion yields are given. They obtained a value of 
0.173 ± 0.017 for the K + /n+ ratio and 0.143 ± 0.014 
for the K~ /ir~ . At midrapidity the density of protons 
plus antiprotons, dN p+ p/dy is about half the density of 
charged kaons, dN K + +K - /dy, as can be read off from 
Fig. 2 of From these one can conclude that about 
80% of the total charged particles are pions. This value 
can be consistent with our previous estimate that 63% 
of total energy is carried by pions, taking into account 



that - due to their larger mass - kaons and protons carry 
more energy on average then pions. 

We also have to check the validity of the increasing 
entropy condition, Eq. (|42p . For that purpose we have 
to specify the value of the Stefan-Boltzmann constant in 
the QGP EOS e = cxT 4 . Assuming a mixture of massless 
quarks/antiquarks and gluons with Nf quark flavors and 
N c colors its value is 

v=-(-J^ + Nt-l). (53) 

Using Eq. ([3]) we get for the entropy density in the QGP 

dP 4 . 
sqgp =gf= g^ 3 - (54) 

We assume a quark-gluon plasma with two quark flavors 
N f = 2, and iV c = 3. 

Because the four-velocity is unchanged during the 
phase transition, the entropy condition (I42[) is equivalent 
to [s] > 0, or spfo/sqgp > L where sqgp and spfo de- 
note the entropy density valid at the quark-gluon plasma 
and post freeze-out side of the FO hypersurface, respec- 
tively. Figure^] (a) shows the ratio spfo/sqgp as a func- 
tion of the temperature in QGP. In order to obtain this 
quantity one has to calculate the post FO temperature 
in terms of the pre FO (QGP) temperature. This is done 
using the QGP EOS e = cT 4 and inverting numerically 
the PFO EOS Eq. (gU). Then the post FO entropy den- 
sity is given by Eq. (J46J) . 

The continuous curve on Fig. [5] (a) shows that the as- 
sumption of the simple pion gas on the hadronic side of 
the phase transition violates the entropy condition, giv- 
ing spfo/sqgp ~ 0.1 — 0.15. To understand the reason 
of this we recall that the entropy carried by free massless 
particles in an equilibrium ideal gas is a universal con- 
stant (3.6 for bosons, 4.2 for fermions), thus, the entropy 
of the system is proportional to the number of parti- 
cles. For massive particles the entropy per particle in- 
creases but this increase is not very large for pions (see 
Q). On the other hand the number of particles decreases 
significantly. First of all because the degeneracy of states 
is reduced drastically during the phase transition (from 
AN F N C = 24 for quarks plus Nq - 1 = 8 for gluons 
to 3 for pions). Therefore, the available energy has to 
be distributed among higher energy states allowing for a 
smaller number of created particles. (Also, some of the 
available energy is needed to generate the pion mass.) 

On Fig. [2] (a) we also plot the entropy density ratio 
for two other possible models of the hadronic medium 
already mentioned at the end of Section IIVI Including 
additional hadrons (curve "hadron gas" on the plot) we 
increase the degeneracy of the states, and also, higher 
mass hadrons carry more entropy, but the increasing en- 
tropy condition is still not fulfilled. On the other hand, 
the constituent quark/ antiquark + pion gas model is con- 
sistent with the increasing entropy condition for QGP 
temperatures > lAQMeV, which is fulfilled for the part 
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of the FO hypersurface where y < 3, as can be seen on 
Fig.[2](b). At the edges, where the fluid rapidity is larger, 
FO happens too late in the Landau model. The increas- 
ing entropy condition favors an earlier FO at the sides of 
the system. 



On the other hand, the constituent quark/antiquark 
+ pion gas model describes an intermediate state of the 
system where hadrons other than pions have not yet been 
formed. Exactly, the coalescence of quarks to hadrons is 
the process responsible for the reduction of the number 
of particles, and, thus, for the usual problems with the 
increasing entropy condition. 
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y 



FIG. 2: (a) (Color online) The ratio of entropy densities in 
the hadronic (or post freeze-out; PFO) state and the QGP 
state as a function of the temperature in QGP. The solid line 
represents the results obtained assuming a simple pion gas 
on the hadronic side of the phase transition. The other two 
curves show the results of the "hadron gas" and "constituent 
quark/antiquark + pion" gas models, (b) Temperature in the 
QGP phase at FO in the Landau model as a function of fluid 
rapidity. Parameters of the initial state correspond to Au+Au 
collisions at ^/snn = 200 GeV, as described in the text. 



Carrying out the calculation of the pion rapidity dis- 
tribution assuming the hadron gas model on the post FO 
side one finds that the shape of the pion rapidity dis- 
tribution is distorted: there is an enhancement of pion 
production at large rapidities. The reason for this non- 
realistic pion rapidity distribution is lying in the assump- 
tion of thermal and chemical equilibrium in the post FO 
medium and the r = const, choice of the FO hypersur- 
face. The FO temperature, that determines the ratio 
of produced particles of different mass, is not constant 
along the FO hypersurface (see Fig. [2] (b)). At the side 
of the fluid with large fluid rapidity the temperature is 
low, therefore a smaller number of large mass hadrons 
is created and more energy is left for pion production. 
This explains the enhanced pion production at large ra- 
pidity in this model. A key element in this effect is the 
presence of large differences between the masses of the 
particles present in the post FO medium. 



We can conclude that the present hadronization de- 
scription in the Landau model can give a good descrip- 
tion of the pion rapidity shape only if the particles in the 
post FO medium have similar masses. 
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